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ABSTRACT 

IChang et al.l (j2006l ) reported millisecond duration dips in the X-ray intensity of Sco X-1 and attributed them to 
occultations of the source by small trans-Neptunian objects (TNOs). We have found multiple lines of evidence that 
these dips are not astronomical in origin, but rather the result of high-energy charged particle events in the RXTE 
PCA detectors. Our analysis of the RXTE data indicates that at most 10% of the observed dips in Sco X-1 could be 
due to occultations by TNOs, and, furthermore, we find no positive or supporting evidence for any of them being due 
to TNOs. We therefore believe that it is a mistake to conclude that any TNOs have been detected via occultation of 
Sco X-1. 

Subject headings: X-rays: general. X-rays: individual (Sco X-1), solar system: general. Kuiper Belt 



1. INTRODUCTION 

IChang et al.l ()2006l ) found statistically significant 1-2 
millisecond duration dips in the count rate during X- 
ray observations of the bright X-ray source Sco X-1 car- 

[ ried out with the Proportional Counter Array (PCA) on 
the Rossi X-ray Timing Explorer (RXTE) and attributed 
them to occultations of the source by small objects or- 

[ biting the Sun beyond the orbit of Neptune, i.e., trans- 

'. Neptunian objects (TNOs). In all. IChang et all ()2006l ) 
found some 58 dips in approximately 322 ks of Sco X- 

[ 1 observations. Given that the RXTE spacecraft moves 
through the diffraction-widened shadows of any TNOs at 
a velocity of ~ 30 km s~^, dips of ^2 ms duration should 

I correspond to a TNO size of ~60 m. If the identification 
of these dips with occultations by TNOs is correct, the 
dips would provide extremely valuable information on 

\ the number and distribution of solar system objects of 
~ 20-100 m in size. We have found evidence that these 
dips are produced by electronic dead-time as a result of 

] high-energy charged particle events in the RXTE PCA 
dete ctors. Prel i minary repo rts of our results were given 
by J ones et al.1 (|2006l l2007f ): herein we give a more de- 

] tailed and complete report. 

2. AVERAGE PROPERTIES OF DIPS IN THE SCO X-1 
COUNT RATES 

Subsequent to the report bv IChang et al.l ()2006f ). we 
searched for dips of the type they describe in the archival 
data obtained in ~ 880 ks of RXTE/PCA observations 
of Sco X-1. These observations were performed starting 
in early 1996 soon after the launch of the spacecraft. 
For the search we only used observations that provided 
count rates of events with a time resolution of < 0.25 ms. 
The total count rate depends on the strength of Sco X-1 
at the time of observation as well as on the number of 
Proportional Counter Units (PCUs) that were operating 
and the location of the source within the field of view; it 
ranges from below 40,000 cts s~^ to nearly 200,000 cts 
s^^ (see Figure [l]). 

Coincidences within a ^ 10 fis window among two 
or more of the measurement chains in a Proportional 
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Fig. 1. — Cumulative distribution of exposure time as a function 
of PCA count rate during observations of Sco X-1 which yielded 
high time resolution data and were used for the present analysis. 
The ordinate gives the total exposure time at count rates greater 
than that of the corresponding rate on the abscissa. Two-LLD 
events have been included in the count rates when available (see 
text). 

Counter Unit (PCU) are normally used to identify 
charged particle events (see lJahoda et al.ll2006l . and ref- 
erences therein, for technical information on the PCA). 
However, the intensity of Sco X-1 is so high that there 
is a substantial count rate due to the detection of two 
X-ray photons in two distinct regions of the detector ser- 
viced by different measurement chains within the 10 //s 
window. For most of the Sco X-1 observations, the rates 
of such so-called two lower- level discriminator ( "2-LLD" ) 
events were telemetered with sub-millisecond time reso- 
lution. In those cases where the 2-LLD event data are 
available, we add 2 counts for each 2-LLD event to the 
counts of single LLD events. 

In our search, we looked for instances when the number 
of counts in a time "window" that is an integral multi- 
ple of 0.25 ms (2^^^ s to be precise) was less than that 
reasonably expected to occur given the mean count rate. 
The search was performed for each of the seven window 
intervals of 0.75, 1, 1.5, 2, 2.5, 3, and 4 ms. The expected 
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number of counts in each window was determined from 
the running average count rate in a time interval cen- 
tered on the time window of 128 ms for the 0.75 and 1 
ms windows, 192 ms for the 1.5 ms window, or 256 ms 
for the 2, 2.5, 3, and 4 ms windows. 

Given TV, the number of counts actually measured in 
the window, and N^xp, the expected number of counts 
in that window based on the running average, we com- 
puted P{n < N\Nf.xp), the probability that N or fewer 
counts would be detected based on the simplistic assump- 
tion that the counts obey a Gaussian distribution with 
mean Ng^p and standard deviation ^JN^. We define 
the detection of a dip as any instance in which we found 
P(n < N\Nexp) < 10^^°. In the search, the same dip 
could be found multiple times, e.g., in contiguous inter- 
vals when it was longer than the window interval or for 
different window intervals. We generated a list of 203 
dips in which these duplications had been eliminated. 
Of these, 196 occur in data in which 2-LLD events had 
been taken into ac count. All but three of the 58 dips of 
IChang et all (l2006f) and aU but 14 of the 107 dips identi- 
fied bv lChang et all (l2007l ) were identified in our search. 
Eight of the latter 14 dips were not fou nd by us be- 
cause we did not search 5 orbits of data that lChang et al.l 
()2007t ) searched, and because 6 of the 14 did not meet 
our significance threshold (perhaps because we included 
2-LLI ) events). We detected more dips than lChang et al.l 
()2007ft because we included 2-LLD events, because our 
significance threshold was slightly lower, and possibly be- 
cause of other minor differences in the searches. 

The frequency of occurrence of dips was found to be 
a function of count rate, and, as one should expect, 
tended to be higher at the higher count rates. Given 
our detection criteria, we found no dips at count rates 
below 43,000 cts s""'^ and only three at count rates be- 
low 55,000 cts s~^. We searched for dips in data which 
included approximately 570 ks of observations in which 
the count rates were above 55,000 cts s~^ (Fig. [1]). For 
each of the seven window durations, we computed N and 
Nexp for window intervals that started every 0.25 ms. 
Thus the number of independent trials must be less than 
7 X 570,000/0.00025 - 1.6 x 10^°. If the actual proba- 
bility of finding a dip due to a statistical fluctuation in 
each trial is < 10"^", then at most a few of the 203 dips 
could be the result of statistical fluctuations. 

In every case, the actual probability of getting or 
fewer counts differs from our computed value for at least 
two reasons. First, the Gaussian distribution we uti- 
lized overestimates the probabilities of small numbers of 
counts by large factors when compared with a Poisson 
distribution of the same mean. Second, the intensity of 
Sco X-1 is time variable, and there is some chance that 
the mean intensity at the time of a dip was lower than 
that corresponding to N^xp- If source variability was 
important, our detection procedure would have led to 
numerous spurious detections, particularly for the longer 
window durations. However, only four dip-like events 
were most significantly detected when using either the 3 
or 4 ms window durations. We further checked the ef- 
fects of source variability by, first, integrating properly 
normalized power density spectra of Sco X-1 count rate 
data over the frequency range of 5-300 Hz, 4-250 Hz, 
4-200 Hz, 4-170 Hz, or 4-125 Hz for the 1.5, 2, 2.5, 3^ 
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Fig. 2. — Distribution of fitted RMS widths and minimum nor- 
malized count rates for 203 detected dip events. The minimum 
normalized count rate is the value of the model count rate at the 
center of the fitted Gaussian divided by the average count rate 
away from the dip. 

and 4 ms window durations, respectively. The resulting 
root-mean-square fractional amplitudes in the given fre- 
quency bands are typically ~ 2.5% and are almost always 
less than ~ 4.5%. These variablity estimates for each of 
the five window durations were then used together with 
the average count rates on an (RXTE) orbit-by-orbit ba- 
sis to estimate the effects on the dip detection probabil- 
ities. The results are consistent with the paucity of dips 
detected most significantly using the 3 ms and 4 ms win- 
dow durations and indicate that statistical fluctuations 
in the presence of source variability did not produce more 
than ^ 5 spurious detections. 

In order to obtain estimates of dip widths and depths, 
we used the IDL procedure gaussfit to fit each dip profile 
with a function f(t) representing a constant count rate 
plus a Gaussian shaped dip: 

/(f) = A-Se-(*-*°)'/2<T^ (1) 

where A, B, Iq, and a were the parameters to be deter- 
mined. The values of the fitted widths as parameterized 
by the values of ct, and of the fitted minimum normalized 
count rates, i.e., the values of 1 — B/A, of the ~ 200 dips 
are shown in Figure [2l The key features of this figure 
are the relatively narrow range of widths and absence of 
long dips. Nearly all of the events have 0.4 < a < 0.8 
ms, and there is only one dip with cr > 1.1 ms. These 
aspects of the plot are discussed below (see Section |4]). 

If the counting rate dips arc the result of occultations, 
then we would expect that diffraction effects would pro- 
duce small count rate increases on either side of the dips. 
The sizes of these sidelobes and indeed the other details 
of the dip profiles depend on a number of factors in- 
cluding the sizes, shapes, distances, and velocities of the 
occulting bodies, the impact parameters characterizing 
the occultation events, and the velocities of the RXTE 
spacecraft at the times the dips were observed. 

To quantify our expectations of the average diffraction 
sidelobe size and other typical properties of dips pro- 
duced by occultations by TNOs, we performed a Monte 
Carlo computer simulation of an ensemble of occultation 



3 




Fig. 3. — Average profile of the 2189 detected 25 km/s model 
population occultation events. The profiles are superposed after 
centering and stretching according to the parameters of the fitted 
Gaussian curves and have been normalized to the count rates at 
more than 0.007 s from the dip center. 



events in which the occuhing bodies were assumed to be 
opaque spheres, to have radii s that foUow a distribution 
^ ~ s""*, and to foUow prograde circular orbits at a dis- 
tance of 40 AU from the Sun. Since the effects of diffrac- 
tion are wavelength dependent, we took the spectrum of 
detected X-rays to be that of a typical pulse height spec- 
trum of Sco X-1 measured with the PCA. The simulation 
was carried out for each of four cases in which the rela- 
tive velocity between the spacecraft and the shadows of 
the occulting bodies was a fixed value, viz. 15, 25, 30, or 
35 km s~^. 

Simulated dip profiles were computed for a large num- 
ber of occultations with object size and occultation im- 
pact parameter chosen at random according to the ap- 
propriate distributions. The profiles were normalized to 
unity at times far from the occultation centers and were 
then inserted into the real PCA count rate data by mul- 
tiplying the actual count rates by the dip profiles. The 
dip center times were chosen at random among the 880 
ks of observations that were used for the present analysis. 
These data were then searched for dips using the search 
algorithm described above. 

We fit the detected model dips with the function given 
by eq. (1). The light curves of the PCA data containing 
the detected simulated dips were superposed after align- 
ing the dip centroids and rescaling the time scale of each 
light curve in order to normalize all of the fitted widths 
to the value a — 0.85 ms. Interpolation of the rescaled 
bin times to 0.25 ms time bins was required to superpose 
the rescaled light curves, and, as a consequence, adjacent 
bins in the superposition are not completely statistically 
independent. The results are shown in Figure [3l 

The average profile of the real detected dips was simi- 
larly constructed by superposing the PCA light curves of 
202 of the 203 detected dips after alignment and stretch- 
ing; the one dip with width a ^ 2.3 ms was excluded. 
The results are shown in Figure ID To check this, we 
also superposed only those 109 dips with fitted widths 
in a narrow range, i.e., with 0.55 ms < a < 0.75 ms. 



Fig. 4. — Average dip profile of 202 dip events, normalized to the 
mean count rate at more than 0.007 s from the dip center (dashed 
line). The standard deviation computed for these bins is 0.016 
(1.6% of the count rate). 

The light curves containing those 109 were aligned and 
summed, but no stretching was done. The results are 
shown in Figure [5] 

Per the results shown in Fig. [31 we expected to see 
sidelobes with intensities as high as ~8% above the mean 
count rate determined substantially away from the super- 
posed dips. The superposed light curves do not exhibit 
diffraction sidelobes as high as those evident in the av- 
erage profile of the simulated dips, despite having statis- 
tics sufficient to reduce fluctuations to ^1.6% (1 a) of 
the mean count rate. While the differences are not suffi- 
ciently significant to be conclusive, they strongly suggest 
there may be a problem with the occultation hypothesis. 
We elaborate on this in the discussion section below. 

Three additional potential problems with the occul- 
tation interpretation are manifest from the dip profiles. 
First, the summed dip profile is distinctly asymmetric in 
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Fig. 5. — Superposition of all 109 events with width parameters 
between 0.55 and 0.75 ms, with 0.25 ms time bins, normalized 
to the mean count rate at more than 0.007 s from the dip center 
(dashed line). 
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shape as IChang et al.l l)2006f ) suggested for many of the 
individual dips. Second, the distribution of dip widths is 
narrower than what one would expect from occultations 
by bodies with a power-law size distribution of index -4, 
i.e., there are fewer than expected statistically significant 
dips with Gaussian FWHM widths greater than ~ 1 ms. 
Third, when diffraction effects are taken into considera- 
tion, one would expect to see a correlation between the 
fitted widths and the fitted minimum count rates, such 
that longer dips tend to be deeper on average. No such 
correlation is seen in Figure [2] 

3. SEARCH FOR AN ALTERNATE EXPLANATION 

These findings prompted us to further explore alter- 
native explanations for the dips. Only one hypothesis 
appeared to be worthwhile to pursue, i.e., that the dips 
are caused by electronic dead time in response to some 
type of charged particle shower in the spacecraft. Un- 
fortunately, no information with millisecond time resolu- 
tion was available on the non- X-ray background during 
the Sco X-1 observations. Counts of good events, very 
large events (VLEs), propane-layer events, and a catch- 
all category of other types of events (hereafter SMI-other 
events) that includes multiple LLD events are available 
at 1/8 s time resolution from Standard Mode f; most 
other types of data are only available with 16-s time res- 
olution. 

The VLB flag for a PCU is set when the electronics 
detect an event in that PCU with energy greater than 
~100 keV; this can happen in response to the ionization 
produced by a single charged particle or to that pro- 
duced by multiple charged particles which penetrate the 
detector nearly simultaneously. Such a large event can 
produce ringing in the front-end of an electronic measure- 
ment chain. Therefore, in response to the occurrence of a 
VLE, the digital logic shuts down the electronic process- 
ing of events in that PCU for a fixed time period, chosen 
to be 50 /iS for almost all of these Sco X-f observations. 
In addition, each of the 6 main xenon-layer measuring 
chains is disabled until its charge drops to an acceptable 
level. In order for the detector to be shut down for an 
extended period (> 1 ms), an extraordinary amount of 
charge must be deposited on most of the 6 main mea- 
suring chains; it is unclear, at present, whether this can 
happen in response to a single charged particle. 

Figure [H] shows counts of three different types of events 
from Standard Mode 1 in 1/8-s time bins superposed 
around the times of 201 dips. In each panel, the centers 
of the dips have been placed in the bin at time = 0. The 
top panel shows the rates of good events, i.e., those not 
identified as being due to charged particles, in the main 
xenon layers of all operating PCUs, and clearly shows 
the superposed dips; two-LLD events are not included 
in these rates. The counting rate drops by only ~0.7% 
because of the dilution of a ^ 2 ms dip within a 128 ms 
bin. In contrast, the middle panel shows the enhance- 
ment of the counting rate of SMI-other events in the 
vicinity of the dips. The peak is highly significant (~38 
a). The bottom panel corresponds to the VLE event 
rate superposed around the dip times. This peak is also 
statistically very significant (~7 a). The increase in the 
VLE rate is nearly so large as to be consistent with the 
detection of ^1 VLE per PCU per dip. We discuss this 
further below. 



The enhancements in the SMI-other event and VLE 
rates around the times of the dips indicate that there is 
an increase in the rate of detection of non- X-ray events. 
We speculate that these non-X-ray events interrupt nor- 
mal event processing for 1-2 milliseconds in most of the 
PCUs roughly once per hour due to the collection of very 
large amounts of charge. Such an energetic event may be 
the consequence of a particle shower produced by the col- 
lision of a high-energy cosmic ray with a nucleus in the 
RXTE spacecraft. In any case, further clarification of 
the causes of the observed dips would be of interest. 

4. DISCUSSION 

The observed dips have widths and depths that are 
approximately what one might expect to be produced 
by occultations by TNOs, even though much wider 
dips would be detectable in principle (given appropri- 
ate depths). Thus we are obligated to seriously consider 
the hypothesis that some or all of the observed dips are 
the product of TNO occultations. However, close exam- 
ination of the RXTE PCA data reveals six signatures 
that independently indicate that few and possibly none 
of the observed dips are due to occultations by TNOs. 
The signatures are (1) the numbers of SMI-other events 
during the dips; (2) the numbers of VLE events during 
the dips; (3) the absence of the expected diffraction side- 
lobes; (4) the temporal asymmetry of the dips; (5) the 
almost total lack of dips longer than ~1 ms; and (6) the 
lack of correlation between dip duration and depth. We 
discuss each of these in turn. 

(1) SMI-Other Events: Fig. [6] shows that there is, 
on average for 201 dips, a large excess of SMI-other 
events at or near the times of the dips. On average, in- 
dividual dips should show an excess of SMI-other events 
at the 2.7(7 level. In Fig. [7] we show a histogram of the 
number of SMl-othcr events in the 1/8-s time bins cor- 
responding to the dips expressed in standard deviations 
from the mean. The mean value obtained from the his- 
togram is ^ 2.7(7, as expected. 

If one makes the reasonable assumption that the num- 
bers of SMI-other events should not be affected by true 
occultations (other than by negligible increases due to 
reductions in the electronic dead time) , then one may es- 
timate the maximum fraction, /occ, of the observed dips 
that represent genuine occultation events that is consis- 
tent with this distribution of SMI-other events. We con- 
structed the following simple function with which to fit 
the histogram, and thereby constrain face- 

probability = ^/_e"^^/2+— L=(l-/,,,)e-(^-^-)V2-^^- 

(2) 

where C is the number of excess SMI-other counts (in 
units of standard deviations of the counts per bin in each 
PCA light curve), and Ccr is the mean of C for those dips 
which are not the results of occultation events and which 
we take to be « 2.7/(1 — /occ)- The distribution of the 
numbers of excess SMI-other counts is wider than what 
would be expected from a Poisson distribution with a 
mean equal to the slightly increased (on average) num- 
ber of SMI-other events per bin; the width of this com- 
ponent of the fitting function is adjusted by means of the 
parameter cr^r. Fits of the function to the histogram in 
Fig. [7] were carried out with fits using both Gaussian 
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Fig. 6. — Counts per 1/8-s time bin of different types of PCA detector events superposed, i.e., averaged, around the times of 201 dips. 
The data used for this figure were recorded in Standard Mode 1. The superposition was accomplished such that the bin at time = s 
includes the identified dips. The counts include events from all (typically 3 to 5) of the operating PCUs. Top panel: good xenon counting 
rate data. The decrease in counting rate due to the dips is apparent. The small ('^0.7%) drop in the counting rate is explained in the text. 
Middle panel: counting rate data of SMI-other events that include multiple-LLD events (see text). A highly significant enhancement in the 
vicinity of the dips is evident. Bottom panel: VLE event rate data superposed around the dip times, also showing a statistically significant 
peak. Note that the peak in the VLE event rate is approximately one VLE event per detector per dip event (i.e., ~4 excess events per dip). 



and Cash (1979) statistics. If wc neglect the tail of the 
distribution at high numbers of excess SMI-other events, 
i.e., at > 9(7, we obtain formally acceptable fits with val- 
ues of /occ in the range 0.0 to 0.12 and values of o'er in 
the range 1.85 to 2.35 (based on Gaussian statistics; the 
limits represent the formal joint 95% confidence range). 
Using Cash statistics, we obtain formally acceptable fits 
with values of /occ in the range 0.0 to 0.11 and values of 
(Tcr- in the range 1.65 to 2.15 (95% confidence). These re- 
sults indicate that fewer than 11% of the 203 dips might 
be the product of TNO occultations. 

(2) VLE events: Figure [6] also shows that there is an 
excess of VLE events around the times of the dips. The 
difference between the background rate and that in the 
1 / 8-s bin containing the dips is very close to 4 (actually 
3.9±0.5) extra VLE events per dip. The peak in Fig.[6]is 
significant at the 7a level. If there is precisely one VLE 
per operating PCU for each non-TNO dip, then we would 
expect on average an excess of 4.67 VLEs per non-TNO 
dip. If only the non-TNO dips contribute to the excess 
VLE events then there is an upper limit to /occ that is 
consistent with the observations. If we further allow that 
the statistical mean excess number of VLE events per dip 



may have been as small as 2.9, then a simple calculation 
gives the limit /occ ^ 0.38 (95% confidence). This limit is 
weaker than for the SMI-other events, and, furthermore, 
is compromised by the possibility that more than one 
VLE event could be produced in a operating PCU in a 
single cosmic-ray induced dip. 

(3) Lack of diffraction sidelobes: In Fig. [3] we 
showed an average profile of simulated dips that had been 
inserted into actual PCA data. It should be compared 
to averages of the actual measured dip profiles in Figs. [3] 
and[5l The average model dip profile shows a clear bump 
of ~ 8% amplitude on either side of the dip due to diffrac- 
tion, whereas the averages of the actual profiles show no 
significant evidence for diffraction sidelobes. Thus, we 
conclude that the fraction /occ of legitimate TNO occul- 
tations can be no larger than ~30%, otherwise diffraction 
sidelobes likely would have been detected. Again, while 
this is a clear strike against the dips being due to TNOs, 
the limiting statistically significant constraint that can 
be set due to the lack of diffraction sidelobes is not as 
significant as for the SMI-other events. 

(4) Asymmetry: A comparison of the simulated with 
the actual dip profiles (as in [3]) above, clearly shows 
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Fig. 7. — Histogram of the numbers of excess SMI-other events in 1/8-s time bins corresponding to the times of 201 of the 203 dips. 
The number of excess SMI-other events is given in terms of the square root of the mean number of SMI-other events per bin away from 
the time of the dip. The solid smooth curve represents the best fit with /occ = 0, (Tcr = 2.1 and the dashed curve represents a formally 
unacceptable fit with /occ = 0.24, Ccr = 1.8 (see text). 



a marked asymmetry for the real dip events. This is 
physically implausible if the dips are the product of oc- 
cultation events and therefore testifies against a TNO 
origin for most of the dip events. We estimate that the 
statistical significance of the asymmetry is '--^ 6 cr. Unfor- 
tunately, there is no direct way to use this information 
to constrain the fraction of legitimate TNO occultations. 
The problem is that we do not know, a priori, how large 
the asymmetry is, on average, for non-TNO dips. There- 
fore, we can not tell how 'diluted' the non-TNO events 
are by potentially real ones. Noncthlcss, this marked 
asymmetry is another solid indication that few of the 
dips arc the product of TNO occultations. 

(5) Lack of dips longer than ~1 ms: From Fig. [2] 
we can sec that all of the dips, except for a single event, 
have RMS widths cr < 1.1 ms. In Section [2] we described 
a computer simulation of the production, detection, and 
analysis of dips caused by TNO occultations. For a rela- 
tive speed between the RXTE satellite and the shadows 
of the putative TNOs of w,.g1 = 25 km s~^ we find that 
the fraction of recovered simulated dips with cr > 1.1 ms 
is ^27%. For v^^i = 35 km s~^, 9% of the dips have 
(T > 1.1 ms. We estimate that the average relative veloc- 
ity between RXTE and the shadows of any TNOs was 
not higher than ~ 30 km s~^. For this speed, 16% 



of the dips are characterized by a > 1.1 ms. Therefore, 
if all of the dips are the result of TNO occultations the 
number of longer-duration dips should be ~30, whereas 
the observed number is actually 1. On the other hand, if 
only 15% of the dips are due to TNO occultations, then 
we would expect only ~ 5 dips with cr > 1.1 ms. This ex- 
pected number is marginally statistically consistent, i.e., 
at ~5% confidence, with the detection of one dip with 
a > 1.1 ms. Therefore, we conclude that the lack of 
longer dips allows an upper limit of 15% to be set on the 
fraction, /occ, of potentially real TNO occultations. 

(6) Lack of correlation between width and 
depth: If the dips were due to TNO occultations of Sco 
X-1, we would expect a strong correlation between the 
widths of the dips and their depths. This results from 
the fact that diffraction produces shallow occultations for 
the smaller size occulters, while it produces deeper more 
geometric-shadowing-like occultations for the larger oc- 
culters. As can be seen from the distribution of dip 
widths vs. depths in Fig. [51 there is no such correla- 
tion, with almost all of the dips confined to a narrow 
range of widths (between 0.4 and 0.8 ms) and depths that 
range all the way from 45% to nearly 100%. Thus, the 
fact that the dips we detect include a significant number, 
i.e., ~20%, that are both narrow (ct < 0.7 ms) and deep 
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(minimum normalized count rate below 0.2) whereas only 
~ 2% of the 'detected' simulated dips (for Wrci ^ 30 km 
s~^) are this narrow and deep, indicates that < 10% of 
the dips might be due to TNO occultations. Given the 
effects of statistical fluctuations on the observed number 
of narrow deep dips and the fact that the simulation is 
based upon somewhat uncertain parameters, it is more 
reasonable to use these numbers to set an upper limit of 
- 20%. 

Summarizing the results from approaches (1) through 
(6) above, we find limits on the fraction of valid TNOs 
to be /occ < 11%, < 38%, < 30%, < Q%, < 15%, < 
20%, respectively, where "Q" denotes that a formal limit 
could not be set, but the approach provides an important 
independent indication that the dips are, for the most 
part, not the result of TNO occultations. 

We believe that the combined upper limit on /occ due 
to the joint application of all six approaches is simply the 
minimum value achieved by the most sensitive of these, 
i.e., the constraints cannot be combined. The reason, in 
short, is that the effects we explore serve only to statis- 
tically limit the number of events which could be due to 
TNOs rather than to identify specific qualifying events. 
Therefore, our final limit is simply /occ 5-, 10% . 

One might argue, as did lChang et al.l (j2007[ ). that since 
~ 10% of the observed dips cannot be formally elimi- 
nated as being due to TNOs, they serve as viable poten- 
tial candidates for TNO detections. However, we argue 
that if 90% of the dips can be securely eliminated as TNO 
occultations, and there are six different and independent 
indicators that point in the direction of a common cause 
due to cosmic ray interactions in the detector, then it 
is most plausible that all of the dips have this common 
origin. 

While our results cast serious doubt on whether any 
true occultation events have been detected, one cannot 
yet conclude with a high degree of confidence that no 
such events have been detected. Further investigations 
of the dip phenomenon and its possible causes would be 
of interest. We are working to obtain a new measurement 
of, or upper limit on, the rate of occurrence of occulta- 
tions of Sco X-1 by analyzing the data that are being 
obtained in a new series of RXTE observations of Sco X- 
1 with high-time-resolution information on VLB events. 

5. IMPLICATIONS FOR THE POPULATION OF TNOS 

Given that we detected 203 dips in data that covered 
570 ks of observations with count rates > 55, 000 cts 
s^^, our upper limit on /occ corresponds to an upper 
limit of < 20 occultations in the 570 ks of observations. 
If we adopt a model of the TNO characteristics, this 
upper hmit of detected TNO occultations can be used 
to establish an upper limit on the abundance of small 
TNOs. For this purpose we use the assumptions and 
results of our model TNO simulations described above 
in Section [2l 

A TNO of radius s appears to sweep out a solid angle 
per unit time 

dn ^ 2(S + 5)Vrol 

dt ~ 



15 



(3) 



where s -t- 5 is the maximum impact parameter for which 
this body may produce a detectable dip, Vi-d is the appar- 
ent transverse velocity of the TNO, and D is the distance 
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Fic;. 8. — Summary of measurements of the cumulative size 
distribution of small TNOs including the upper limit determined 
herein from RXTE PCA observations. The HST and the ground- 
based optical measurements to the right of the H ST point have 
been adapted from Table 2 of lBernstein et all 1)20041 ) by converting 
limiting R magnitudes to radii by simply scaling from the HST 
limiting magnitude which we take to correspond to a radius of 
10 km at a distance of 40 AU. We show error bars that give the 
ranges of values consistent with Poisson st atistics at confidence 
levels of 90%. The measurements reported in lRogues et al.| (120061 ) 
arc shown as an upper limit. The lino represents a differential size 

distribution of cx; s~^. 

as 

from Earth to the TNO. We find from our simulations of 
occultations that, on account of diffraction effects, (5 ~ 7 
m for the model TNOs in the relevant size range. We use 
Uioi ~ 30 km/s and D ^ 40 AU. The solid angle swept 
out in time At is 



ilsn^ ^ —At 

at 



2(3 + S)v,c 



At. 



(4) 



For an ensemble of TNOs of various sizes, the average 
solid angle swept out per TNO of radius s > Smin is 
then 



1 



N{s > Srmn) 

2w,.ciAt 



ilsw(s)-^ ds (5) 
ds 



1)2 



(3sto„i/2 + 6) 



(6) 



where we have assumed that the differential size distri- 
bution of TNOs is given by ^ oc s^^ for s > Smin- In 
our simulations we find Smm ~ 15 m. Using this and the 
above values, one obtains 



.sw ^ - 2.8 X 10"^* sr - 9 X 10"'^ deg . (7) 

Given that as many as 20 TNOs may have been detected, 
the upper limit on the areal density of TNOs is then 

20 119 
N{s > 15 m) < , ~ 2 X 10" deg-^. (8) 

9 X 10^11 deg 

This upper limit may be compared to previous mea- 
surements of the size distribution of TNOs at larger 
radii. Figure [5] summarizes the results of surveys of 
TNOs smaller than 1000 km. The smallest TNOs which 
have been securely detected were foun d in the Hubble 
Space Telescope ACS survey reported bv lBernstein et aTl 
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(|2004f) . in which 3 TNOs of radius s > 10 km were found 
in 0.019 square degrees of sky. We also show the popula- 
tion estimates from g round-based surveys th at are sum- 
mar izcd in Table 2 of iBernstein et all ()2004[) . 



We are grateful to Jean Swank for helpful discussions 
on the technical aspects of the PCA and to Jim Elliot 
for a number of helpful discussions. 



REFERENCES 



Bernstein, G. M., Trilling, D. E., Allen, R. L., Brown, M. E., 

Holman, M., & Malhotra, R. 2004, AJ, 128, 1364 
Cash, W. 1979, ApJ, 228, 939 

Chang, H.-K., King, S.-K., Liang, J.-S., Wu, P.-S., Lin, L. C.-C, 

& Chiu, J.-L. 2006, Nature, 442, 660 
Chang, H.-K., Liang, J.- S., Liu, C.-Y., & King, S.-K. 2007, ArXiv 

Astrophysics o-prints, |arXiv:as tro-ph / 0701850 
Jahoda, K., Markwardt, C B., Radeva, Y., Rots, A. H., Stark, 

M. J., Swank, J. H., Strohmayer, T. E., & Zhang, W. 2006, 

ApJS, 163, 401 



Rappaport, S. 2006, 



Jones, T. A., Levine, A. M., Morgan, E. H. 

The Astronomer's Telegram, 949, 1 
Jones, T. A., Levine, A. M., Morgan, E. H., & Rappaport, S. 

ArXiv Astrophysics e-prints, arXiv:astro-ph/0612129ir2 
Roques, F., et al. 2006, AJ, 132, 8TS 



2007, 



